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An SU (Nf) X SU (Nf) Yang-Mills theory on an extra-dimensional interval is considered, with 
appropriate symmetry-breaking boundary conditions on the IR brane. UV-brane to UV-brane 
correlators at high energies are compared with the OPE of two-point functions of QCD quark 
currents. Condensates correspond to departure from AdS of the (different) metrics felt by vector 
and axial combinations, away from the UV brane. Their effect on hadronic observables is studied: 
the extracted condensates agree with the signs and orders of magnitude expected from QCD. 


I. HOLOGRAPHIC QCD 

Since the pioneering work of Q, some attention has 
been drawn to what is called holographic QCD. This 
was further studied using 5D models in ill 11 Holo¬ 
graphic QCD tries to answer the following question: 
If it exists, how does the 5D dual of QCD look like? 
The astonishing success of this bottom-up approach war¬ 
rants a more detailed study of the qualities the dual 
should sho"^ independentl y o f a full-blown stringy set¬ 
up iillllilSlIBiSlIl 

Our modus operandi consists of getting a handle on the 
known behavior of QCD and learning about the features 
the dual 5D model should have. At low energies, one 
possesses theoretical and experimental information con¬ 
cerning pion interactions. In the perturbative regime, 
the OPE of QCD, combined with the large-fVc expansion 
and the softness of some amplitudes are powerful theoret¬ 
ical guidelines for the 5D model’s asymptotic behavior. 
In the intermediate region, experimental data on spin-1 
resonances, like masses and widths, is available. 

Here, we briefly outline the features of our model. 
Our starting point is the fact that gauge fields living 
in a 5D setup are supposed to have a 4D global sym¬ 
metry dual. On the 4D side, QCD possesses a global 
SU (Nf) X SU {Nf) symmetry, whose spontaneous break¬ 
ing gives rise to pions. The main ingredient of the 5D 
model is therefore the chiral symmetry of QCD promoted 
to a 5D gauge symmetry, as well as its breaking to the 
diagonal subgroup. A priori, the main difference between 
the approach of 0, 0 and that of 0 , is the implemen¬ 
tation of the symmetry-breaking. We discuss later how 
this influences physical predictions. 

In our model we brought into play two ideas related to 
this breaking. First, boundary conditions (BCs) in the 
extra-dimension represent the spontaneous breaking by a 
4D condensate of infinite dimension, before the AdS sin¬ 


gularity has been smoothed out. Second, the gauge sector 
builds up the pion field: the Goldstone boson (GB) na¬ 
ture of the latter is thus ensured (protected to all orders 
by the 5D gauge symmetry |l^). These two inputs were 
enough to provide a massless pion with the right trans¬ 
formation properties: the model can be rewritten as a 4D 
lagrangian whose interactions explicitly obey chiral sym¬ 
metry (for details, see 0). This ensures that Green’s 
function computed in the model will satisfy the Ward 
identities of chiral symmetry. 

This picture of chiral symmetry breaking was just the 
entree: in the massive sector, the two infinite towers of 
vector fields provided by a Kaluza-Klein decomposition 
have just the right symmetries to ensure softness of am¬ 
plitudes at high energies. The reason behind is the under¬ 
lying gauge structure which, from the 4D point of view, 
manifests itself in the form of clever sum rules between 
resonance couplings ii Moreover, these two towers 
provide the partonic logarithm expected from asymptotic 
freedom, as a consequence of the asymptotic conformal 
invariance in the extra-dimension. 

Furthermore, the lowest excitations could be identified 
as the observed p and ai mesons of QGD. Given as in¬ 
puts the pion decay constant and the rho mass Mp, 
the mass of the ai was predicted within the experimen¬ 
tal range. Another conspicuous result of the 5D model 
is that the experimentally tested rho-meson dominance 
is an automatic property, whereas it was an assumption 
in 4D models for resonances. It is again a consequence 
of the extra-dimensional nature of the model, and man¬ 
ifests itself in 4D as the decrease of the couplings of the 
massive resonances to the pion. Since dominance of the p 
is satisfied within a few percent, and the soft high-energy 
behavior is ensured, the 5D models succeed in predicting 
the low energy constants of Chiral Perturbation Theory 
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Such 5D models therefore match well to the low-energy 
predictions of QCD (which depend on chiral symmetry 
and lightest-meson dominance), as well as onto the HE 
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partonic logarithms ^. In this work, we ask in addition to 
reproduce the analytic form of the first terms in the OPE 
as given in large-iVc QCD. Various Ansatze which per¬ 
form this task for two-point functions can be constructed, 
however the 5D approach presents the following appeal¬ 
ing features: 

• There is no need to impose relations between (an 
infinite number of) constants to constrain the HE 
behavior: the 5D gauge symmetry takes care of 
that task. 

• The condensates (power-corrections to the AdS 
metric) are input parameters in the 5D model, 
rather than derived quantities. The consequences 
of such deviations from conformality have com¬ 
putable effects on two-point and higher Green’s 
functions EEia 

We provide the tools for studying the model in Sect ion ITU 
analyze the results in Section inTl and offer the conclusions 
in Sect ion ITVI Appendix IXI proves that our procedure re¬ 
spects covariance and locality and in Appendix^one can 
find the approximate dependence of hadronic observables 
on the condensates. 


II. INCORPORATING CONDENSATES IN THE 
5D MODEL 

The simplest model presented in 0 was useful to il¬ 
lustrate a basic point: imposing the above-mentioned 
symmetries (AdS geometry and 5D gauge SU(V/)i x 
SlJ{Nf)n ^ SlJ{Nf)v) and adjusting two experimental 
inputs (/tt and Mp) was enough to obtain an agreement 
for low-energy quantities at the level expected for a model 
of large-Vc QCD at leading order. 

The next stra is to incorporate more features of QCD: 
the model of 0, defined on truncated AdS space with 
symmetry-breaking BCs lacks power corrections in the 
correlators, known to be present in QCD. Indeed, the 
OPE, for large euclidean momentum = —q^ 0 
should read, assuming factorization in the large-Ac 

nv.^(-Q^) = -^|A + Iog(^)+O(a«)| 

1 as{Gp„Gi^‘') / -I \ 28 nasiqq)'^ 
127r Q4 +1^11/7^ 9 Q6 

where A is a subtraction constant. Note that dimension 
4 and 6 condensates have quite different chiral proper¬ 
ties: the former is chiral-invariant and therefore does not 


contribute to Hy — H^, whereas the latter contributes 
essentially to Hy — H Amt- This is why we will con¬ 
sider both of them. In Q the pieces of the second 

line in ED were absent, and numerical agreement for 
Nc was not imposed. 

From the model building point of view, there are many 
ways to introduce the condensates but, at the end of the 
day, they are just modifications of the vector and axial 
wavefunctions located at a distance from the brane where 
the electroweak sector lives Since the model of 0 al¬ 
ready boasts a pion, there is no need for an extra spin-0 
field: we therefore introduce all condensates as modifica¬ 
tions of the metric. We will then be able to originate both 
axial and vector wavefunction modifications at once in a 
parametrically appealing fashion. This will have further 
consequences in the low-energy physics, since the modifi¬ 
cations will be especially relevant there: we shall examine 
them in Section cni 


A. Modifications of the metric and OPE 


The extra dimension considered here is an interval. 
The two ends of the space are located at Iq (the UV 
brane) and h (the IR brane). With a purely AdS met¬ 
ric w (z) = Iq/z, the space-time interval then reads in 
conformal coordinates 

ds^ = w ( 2 )^ [rjp^dx^dx^ — dz^) . (2.2) 

This reproduces the logarithm piece in JZU, with only 
exponentially decaying corrections. Let us see how the 
power corrections of EU can be generated. The high- 
energy expansion for the two point functions can be ob¬ 
tained by solving the following differential equation for a 
function (fx (z), {X = V,A) 

-dl- (log wx (z))' dz}^x {z) = 0, (2.3) 

where we have used the euclidean 4D momentum = 
—q^ > 0. The boundary conditions to be imposed on ipx 
on the IR brane are the same as for the fields, respectively 
dz^v iQiZ)\z=ii — 0 V’A {Qih) = 0, while the value 
on the UV brane ipx {Q,lo) is just a normalization. The 
two-point function is then obtained through 00111 
0 


TT t ^2^ 2 1 dz^px (Q,z) 

Ax {-Q ) = 2 - 77; — 

ffiQ m (Q, z) 


(2.4) 


z=ln 


With a purely AdS metric, w (z) = l^jz, one obtains for 

Q > l/h 


(f{Q,z) ~ QzKi (Qz), 


(2.5) 


^ Note that this is realized with a fast-growing resonance spectrum: 
the mass of the n-th KK mode goes as n for n ^ 1, rather than 
y/n obtained from linear Regge trajectories. 


2 In iS: couplings of spin-1 fields to a scalar sector were respon¬ 
sible of modifying the axial wavefunction only. No 1/Q® were 
generated in the vector two-point function, in disagreement with 
i2.H . Neither were terms included, also in contradiction 

with 12.li . 
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which reproduces the log piece in with the matching 


N,= 


127r^/o 

5i 


( 2 . 6 ) 


Assume now that the metrics near the UV take the 
form 

Y^wv,a{z) l + ^JO,)iz-lo)^ 

+ o((z-;o)®), (2.7) 


where we have already adopted a convenient name for 
the coefficients {O 4 ) , {Oq). With the metric 12.7II . one 
solves the differential equation to obtain Ilv,A via 

(E3I) 


n^A {-Q^) = 


+ 


Nc 

127r2 






+ 0 




m 

( 2 . 8 ) 


Deviations from conformality with a given power of 
in the metric yield a power of in the two- 

point function, with a computable coefficient. Note that 
0{as) contributions in 12.1 II can be accounted as inverse 
powers of log 2 : in the metric However, we restrict our 
attention to the condensates here. 

In the following, we will use the particular functional 
form 


or smaller. Plugging this back into one sees that 

the deviations of the metric are of order one for the di¬ 
mension 4 and 6 condensates, but decrease quite fast with 
d, as • Comparing our model to data in Section 

nni we will indeed extract values for condensates that 
are within the expected range of 12.1011 and of QCD phe¬ 
nomenological analyses, confirming the consistency of the 
whole approach. 

Note that, since it is an order parameter, the 
dimension-6 condensate generates a power correction to 
Hy — Ha, implying that there are only two vanish¬ 
ing Weinberg sum rules. To reproduce this, it is nec¬ 
essary that the wavefunctions of the vector and axial 
helds feel a different metric. Otherwise one would have 
(Oe)^ = (^6) y in in contradiction with the QCD 

OPE Eg. 12.HI . We explain how this essential point is 
achieved in the next Section and Appendix fXl 


B. Definition of the model 

The model is defined by an SU {Nf) x SU {Nf) YM 
lagrangian in an extra dimension. We limit ourselves to 
terms involving two derivatives. Denoting by wv,a (z) 
respectively the warp factors that will be identified later 
as the metrics felt by the vector and axial fields, we start 
from the action 

Wviz)+WAiz) ^ : T, T, ^ 

- 2 - K^MS^NT + ^MS^NT) 

+ {wv (z) — WA (z)) (LmS^^RnT^)} , (2.11) 



z 

for which a given power of z‘^‘^ in yields ex¬ 
actly a power of only when multiplying 

(E3 by exp ^^ 2 d) {z - is a term 

{ 02 d) /Q^'^ generated in 12.811 . 

Representing deviations from conformality as devia¬ 
tions of the metric from AdS provides a basis to enforce 
the OPE of QCD term by term. Including the lowest 
dimension condensates should be sufficient. Indeed, one 
would naively expect a condensate of dimension 2d to be 
of order the confinment scale to the appropriate power, 
i.e. Aqqq, and to be proportional to Nc in the large-A^c 
limit. Trading Aqcd for the scale l/h in the 5D model 
wich induces the mass gap, one thus expects 

{02d) ~ (2.10) 


^ We thank Ami Katz for pointing out this fact to us. 


where the capital letter index runs over the five dimen¬ 
sions. The last term in (ITnli summarizes the low-energy 
effects of couplings to other fields, and will ultimately be 
responsible for the vector and axial metric being different. 
It involves a unitary auxiliary field Q {x^, z), transform¬ 
ing as a bifundamental under the gauge groups. This 
spurious field is eliminated by imposing the constraint 
Dsfl = 0. 

Note that a term of a similar form, and with simi¬ 
lar effects could be written down in the model of Hi, 
using (instead of the spurion fl) a bulk scalar $. Yet, 
since $ is dimensionful, this term could only appear at 
higher orders, presumably suppressed by a large scale. 
In our subsequent fit to experimental data, we shall how¬ 
ever find that the order of magnitude of this term is in 
agreement with our including it at leading order, i.e. as a 
deformation of the metric rather than as a higher-order 
term. This is directly related to the above estimate of 
Nc/li‘^ as the size of the condensates. Quite the oppo¬ 
site from being a large scale, l/h is the only scale in the 
theory: modifications of the metric are of order one (for 
the lowest dimension condensates). 

Another ingredient in |^, responsible for the presence 
of a well-defined pion in the spectrum, consists of the 
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choice of BCs. At the IR boundary z = li: A/j — = 0 

and = 0. At the UV boundary, classical con¬ 
figurations serve as sources (x) , (x) for the analogue 

of the QCD quark currents, i.e. (x, z = Iq) = (x) 

and R^ {x,z = Iq) = (x). Since the spurion satishes 

D^Q. = 0, we need only one boundary condition, and we 
choose D^Vl{x,z = h) = 0. As shown by a covariant 
treatment (Appendix 0, the action (ICTTi is equivalent 
to writing a lagrangian in terms of vector and axial field 
strengths, but with different metrics, respectively wy and 
WA- 

C. Sum rules 

In 0 , four sum rules involving 4D resonances couplings 
emerged as a consequence of the underlying 5D gauge 
structure. Three of them, together with the relations 

Ls = -3L2 = -6Ai, (2.12) 

ensured the unsubtracted dispersion relation for the vec¬ 
tor form factor, as well as the once-subtracted dispersion 
relation for GB forward elastic scattering amplitudes. 
These sum rules still hold with the present modifications 
of the metric, the same consequences for the behavior of 
the amplitudes therefore follow. 

The only sum rule that does not hold anymore when 
WA ^ wv is the one related to the KSFR II ratio: we de¬ 
rive /^/ (En 9v„MvJ > 3 as long as wy (z)-^^) < 0 
(as expected in order to reproduce the result of that 
Ily — 11^4 < 0). Since this sum rule was the only one not 
related to any high-energy constraint, this has no theo¬ 
retical consequences. From the phenomenological point 
of view, it is known that a KSFR II ratio of two is favored, 
but in a model, this depends on the predicted values of /tt 
and Mp-. the agreement is best judged on the observable 
Fp^TTTi-, see Tabled 

On the other hand, the case wa wy gives rise to a 
new sum rule (derived using the completeness relations 
for the wave-functions of the axial KK modes) 

OO 

(/A„+2V2aA„) = 4(L9 + Lio),(2.13) 

n—1 

in the notation of ^ . Each term in the left-hand side van¬ 
ishes separately when wy = wa- Due to this sum rule, 
the axial form factor Ga {q^) satisfies an unsubtracted 
dispersion relation 

OO j^2 

GA{q^) = Y.fA„{fA„+2V2aA„)j^r^-^^) 

n—1 ^ 

The modification wy ^ wa implies immediately Lq + 
Lio 7 ^ 0, and in general + 2 \^aA„ 7 ^ 0 (at least for 
some n). In addition, since the width is propor¬ 
tional to (/ai + 2 \f2aA^ this is accompanied 


in general by a non-zero decay of oi into Try. For a dis¬ 
cussion of the numerical values of these two quantities, 
see Appendix IbI 

D. Quark masses 

The electroweak interactions live on the UV brane, 
which we will therefore also call the EW brane. These 
interactions induce masses for the quarks via the Higgs 
mechanism of the SM. This in turn breaks the chiral sym¬ 
metry explicitly: since we are only modelling QCD, and 
not electroweak symmetry breaking, we simply consider 
hard quark masses. In the present model as in yPT 
such quark masses can be introduced via a spurious field 
X (a:) transforming as a bifundamental under the chiral 
SU (iVy) X SU (iVy) symmetry of our 4D world. In our 
model, this translates as a 4D bifundamental under a 
subset of the 5D SU {Nf)x SU (iV/), namely the one act¬ 
ing on the EW brane. With such an object, one can write 
down the following term on the UV boundary 

^ (xiA^HxJo)) (2.15) 

where one has to keep in mind that x should be counted 
as equivalent to an object containing two derivatives, and 
must be set equal to the quark mass matrix times a con¬ 
stant Bq (of dimension one). 

The key point is then to use the consequences of the 
constraint = 0 and the rewritings of Appendix \B 
which lead to Q (xJo) = U (x), i.e. the GB matrix. We 
thus recover the yPT term /^/4(x{7^). From this it 
follows, strictly as in yPT, that Bq gives a measure of 
the quark condensate at tree level, according to M 

(qq) = -flBo- (2.16) 

As in yPT again, the term 12.1511 also induces a mass for 
the GBs, which become pseudo-Goldstone bosons. Dis¬ 
regarding isospin breaking effects, this reads 

Ml = Bo{mu + md). (2-17) 

We will use the value = 135 MeV in the analysis of 
Section inii 

To summarize the situation in the present model, we 
see that explicit symmetry breaking by quark masses, 
which is given by a 4D lagrangian on the EW brane, 
works exactly as in 4D xPT, provided one realizes the 
equivalence between ft (x, Iq) and U (x). This is also true 
of the higher orders. 


III. NUMERICAL RESULTS WITH 
CONDENSATES 

Once the model has been rewritten in terms of vec¬ 
tor and axial combinations as in Appendix ID one can 
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proceed following the method of 0| to compute masses, 
couplings and decay widths of the pion, rho and oi. In 
addition to the following couplings can be 

easily tested: the coupling of the p to the vector current is 
tested in Fp^ee, its coupling to two pions in Fp^^^r, and 
the couplings of the ai to the axial current in 
which we compare to the experimental ■ We thus 

fit data by minimizing the RMS error on the six quanti¬ 
ties .u, M„, M„,,, with respect 

to their central values provided in |2^. 

We performed a numerical analysis of these quantities 
in terms of the condensates and find that the favored re¬ 
gion in parameter space is consistent with the estimate 
\{02d) \ ^ One can then in turn predict the val¬ 

ues of the low-energy constants Li appearing in the low- 
energy expansion of the model. 

In the favored region, the approximate expression at 
linear order in the condensates are accurate within a 
few percents for /,r,Mp, Mq^ and the Li constants. We 
present these expressions in Appendix 0 


A. Strategy and results 

In Ref.Q, we followed an approach centered on low- 
energies: we imposed = 87MeV (in the chiral limit) 
and Mp = 776 MeV. This led to agreement within the 
expected range (1/W corrections) for low energy quanti¬ 
ties: the mass of the Oi, decays of the p, and allowed to 
predict the low-energy constants of yPT. On the other 
hand, this exhausted the free parameters of the model, 
and yielded a mismatch in the OPE: the coefficient of the 
logarithm led to Nc = 4.3 in Eq. for the chiral limit. 

Another way to see this is that a pure AdS model with 
Nc = 3 but without condensates naturally underpredicts 
Jtt and the yPT low-energy constants (Table P) . This 
mismatch is the evidence that, in order to interpolate 
between very high and very low energies with better pre¬ 
cision, another ingredient is needed: the condensates. In 
the present paper, we start by imposing Nc = 3, as in 
[3,l3- III the pure AdS case, this leaves us with one free 
parameter li, which we adjust to minimize the RMS er¬ 
ror on the six observables ^J^(50/OY/n, where n is 
the number of observables minus parameters (n = 5) . 
The best fit is obtained for li ~ 3.1 GeV“^, yielding an 
RMS error of 12.5%. This is our input Set A, which 
constitutes a benchmark result without condensates: the 
outputs can be read in Table P 


This is because the decays mentioned involve the squares of the 
couplings appearing in the lagrangian: to treat the couplings 
we are interested in on the same footing as Mp, Mai (which 
appear squared in the lagrangian), we consider the fourth-root of 
the decays NV-^aiv-r ■ It is on these six quantities 

that we expect the precision to be of the same order. Note that 
this procedure is very close to that of P, except that we fit 
directly observables rather than couplings. 


1.2 


1 



0 0.2 0.4 0.6 0.8 1 


FIG. 1: Sketch of the favored parameter space: the x-axis rep¬ 
resents {Oi) in units of 10“® GeV^ and the y-axis represents 
{Oe) in units of 10“® GeV®. 


The next step is to include the condensates. Dimen¬ 
sional analysis showed that only the dimension 4 and 
6 condensates are relevant. The dimension-4 conden¬ 
sate {Oi) accounts for one parameter since it is chiral- 
invariant {Oi) = {Oi)y = {Oi )and the dimension-6 
condensates are chosen to satisfy the relation {Oq)a — 
— ^ {Oe)y, derived from factorization: this adds one 
more parameter. Figure ^ indicates approximately the 
favored region for {Oi) and {Oq), such that the RMS er¬ 
ror can be below 9% for an appropriate value of /i. For 
this case, n = 3. 

Inside this favored region, we picked the point denoted 
by B on the hgure. The outputs for this Set B are given 
in Table n 

In TableP we indicate the result for low-energy observ¬ 
ables, without condensate (Set A) and with condensates 
(Set B). We remind the reader of the correspondance be¬ 
tween our parameters and QCD condensates. We present 
the values of the six low-energy observables taken into ac¬ 
count in the fit, as compared to the values extracted from 
1^. We also give the corresponding predictions for the 
low-energy constants of yPT, as compared to the O (p^ 
fit of renormalized at the scale p = 770 MeV. Note 
that the presence of positive {Oi) and {Og) essentially 
improves the prediction for and the Li constants, as 


® At the level we are working, it does not make sense to consider 
O (p®) fits. 
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TABLE I: Results for the two cases: without and with condensates. The signs = X means that the corresponding input 
parameter is fixed to the value X. 



Set A 

Set B 

QCD 

Unit 

h 

= 3.1 

= 3.1 

~ 1/Aqcd 

GeV® 


= 0 

= 0.5 X 10“® 

7^03 (GG) > 0 

GeV'^ 


= 0 

= 0.7 X 10“® 

^7ra„ {gqf > 0 

GeV® 

u 

72.6 

84.3 

92.4 ± 0.3 

MeV 

Mp 

776 

796 

775.8 ±0.5 

MeV 

Ma, 

1.24 

1.34 

1.230 ± 40 

MeV 

r 

p—*ee 

5.63 

5.80 

7.02 ±0.11 

KeV 

r 

-L p— *7r7V 

160 

124 

150.3 ± 1.6 

MeV 

r 

A T —> a 1 i/t 

0.563 

0.580 

0.414 

meV 

lO^Li 

0.36 

0.45 

0.4 ±0.3 


10^1/2 

0.73 (= 2Li) 

0.89 (= 2Li) 

1.35 ±0.3 


10 ® L 3 

- 2.2 (= - 6 L 1 ) 

-2.7 (= - 6 L 1 ) 

-3.5 ± 1.1 


10 ® Lg 

4.7 

5.5 

6.9 ±0.7 


10 ®I/io 

1 

1 

CD 

-5.4 

-5.5 ±0.7 



is readily seen from the approximate expressions of Ap¬ 
pendix 0 In summary, condensates improve the match¬ 
ing between very low energies (/,r and Li's of yPT) on 
one hand, and very high energies {Nc = 3) on the other 
hand. 


B. Interpretation 

In such a model of large-A"c QCD treated at leading 
order, we expect corrections in \/Nc to parameters in the 
lagrangian. To facilitate comparison with |^, we have 
considered square root of such parameters here. Still, 
experience shows that large-A^c models usually fare better 
than this rough estimate, except for some observables 
in which disagreement is expected (an example here is 
see Appendix IBl where the related case of Lg + 
Lio is also discussed). 

Our benchmark Set A, with vanishing condensates, al¬ 
ready has an RMS error of 12.5%. Still, our analysis in¬ 
dicates a preference for non-vanishing condensates: Set 
B has an RMS error of 8.5 % . Fixing the best fit involves 
just low-energy quantities such as decays and masses, but 
from this set of observables one can infer information 
about the high energy. The output is that the preferred 
signs for the condensates agree with those predicted by 
QCD: {Oi) > 0 and {Oq) > 0. In fact, the absolute value 
of {Oi) is consistent with the values extracted from QCD 
phenomenology (see for instance mEnn , and refer¬ 
ences therein). 

As for the dimension -6 condensate, we can estimate 
its size using factorization. However, since we have 
not explicitly introduced the quark masses and the 
strong coupling constant into the model, as (qq) can 
only be obtained from iriq (qq) using input from QCD. 
One starts from the Gell-Mann-Oakes-Renner relation 


rFns. i2.1 612.1 7H . giving a value of nig (qq). Using knowl¬ 
edge of the quark masses and including corrections, [s^ 
finds (gq)(2GeV) = (0.267 ± 0.016 GeV)^. This yields 
28/97ras {qq)'^ ~ 10“^ GeV®, with quite large errors that 
encompass our favored values of (Oq) (Table ^ and Fig¬ 
ure 0. Note that phenomen ologial QCD studies (see 
for instance 1361.1371 1^ and references 

therein) find in general a larger value for the dimension -6 
condensate, more difficult to reconcile with factorization. 

Finally, note that the favored region we have found 
corresponds to an approximate cancellation between the 
effects of {O4) and (Oq) on wy, while they add up in 
the axial channel. This supports the succes of the model 
of [10 , where condensates were introduced in the axial 
channel only, via a bulk scalar field. 


IV. CONCLUSIONS 

We have started from the simplest model of holo¬ 
graphic QCD and described a way to implement term 
by term the OPE of QCD. The 5D model we start with 
is the simplest in the sense that it is defined in terms of 
Yang-Mills fields only, with appropriate BCs that imple¬ 
ment the spontaneous symmetry-breaking. This model 
has only two free parameters: a scale given by the inverse 
of the extension of the extra dimension, and a parameter 
to be identified with the number of colors of QCD. 

In AdS, distance in the extra dimension measures the 
4D momentum: each condensate, giving an giving a 
power-correction in in the deep euclidean, corre¬ 

sponds in the 5D picture to a power-correction to 
the AdS metric, where z measures the separation from 
the UV brane. We use a formalism involving a spurion, 
equivalent to considering two different metrics for the 
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vector and axial fields. This describes order parameters 
in the QCD language. 

The expansion in modifications of the metric is under 
control. First, we find that the naive estimate for its 
size is indeed verified by data and secondly, the effect of 
these modifications become irrelevant as the dimension of 
the condensate increases. This justifies a posteriori our 
procedure of including deformations of the metric rather 
than higher-order operators. 

By introducing the modifications of the metric, we 
show that the infinite number of Weinberg sum rules are 
reduced to the desired number, i.e. two. Also note that 
the high-energy properties of the amplitudes studied in 
0 are not lost, and in fact, a similar property for the 
axial form factor can be demonstrated. 

In terms of input parameters, this model with con¬ 
densates adds two: the dimension 4 and 6 conden¬ 
sates (C> 4 )v = {Oa)^ and {Oq)v = -7/11 (C> 6 )a- 
We performed a complete numerical analysis of the 
effect of those condensates on hadronic observables: 
/^,Mp,Mai,rp^ee,rp^^^ and We find that 

the presence of the condensates improves the agree¬ 
ment from 12.5 % to 8.5 % , even though the num¬ 
ber of parameters has increased froml to 3. More¬ 
over, the favored region in the parameter space corre¬ 
sponds to non-zero values for both dimension 4 and 6 
condensates. The specific value for (O 4 ), correspond¬ 
ing to l/(127r)as is of ( 0.2 0 . 8 ) x GeV^ 

For (Ge), corresponding to 28/97ras the preferred 

range is (0.4 -^1.0) x 10“^ GeV®, in agreement with the 
signs and orders of magnitude expected from QGD and 
factorization for the dimension -6 condensate. 

We also extracted the corresponding low-energy con¬ 
stants: condensates improve the agreement. In partic¬ 
ular, both Lg -|- Lio and become non-zero, in 

connection with the presence of a dimension -6 conden¬ 
sate. However, if factorization is not largely violated, 
both quantities remain smaller than the experimental 
value. 

Apart from the already mentioned differences between 
our model and the one by we mention a more the¬ 

oretical one: in the case of I 3 , the massless pion is an 
admixture containing a bulk scalar, not protected by the 
gauge symmetry. In their construction, the limit of van¬ 
ishing condensates implies the restoration of chiral sym¬ 
metry in the spectrum: equal masses for the vector and 
axial resonances, and the absence of a massless pseudo¬ 
scalar. In our case, vanishing of local order parameters 
does not imply vanishing of non-local order parameters, 
and therefore chiral symmetry is still broken (the axial 
masses are distinct from the vector ones and the Gold- 
stone bosons remain in the spectrum, with ^ 0, see 
Appendix^. The only way to reinstate chiral symmetry 
in the spectrum is then to send to zero. In practice, 
this requires h 00 , i.e. the vanishing of the mass gap: 
a continuum of states is recovered. The converse of this 
statement, namely that confinment implies chiral sym¬ 
metry breaking, is therefore built into our model. 
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APPENDIX A: SYMMETRY-BREAKING, 
COVARIANCE AND LOCALITY 

In the action (ICTTl . we introduced a term 
{Lmn It is clear that such a crossed term is 

going to affect the vector and axial channels differently, 
and thus introduce symmetry-breaking effects. The role 
of the auxiliary 5D field O is the following: one cannot 
simply write down (LmnRst), because this term does 
not respect the 5D SU (Nf) x SU {Nf) group. To make 
it invariant, it appears that we need Wilson lines going 
from the point 2 ; to /i, but this would violate locality in 
the fifth dimension. This would be forgetting one thing: 
to make the whole procedure covariant, we should in ad¬ 
dition reinstate the coset elements on the IR brane, which 
belong to the group [SU (Nf) x SU (Nf)] / SU (Nf) act¬ 
ing a,t z = h- The question is then the following: is there 
a way to write the appropriate combination of these ele¬ 
ments as a local object? 

The answer to the above turns out to be yes, as we 
show here. For simplicity, we start from the solution to 
the problem, and derive the result: we introduce an aux¬ 
iliary field n (x, z) in the bulk, without kinetic term. Pro¬ 
vided the right constraints are applied, this object can be 
decomposed as the above-mentioned product of Wilson 
lines and IR-brane coset elements. Therefore, a spurion 
U with a constraint solves the difficulty mentioned above, 
by allowing a writing that is both covariant and local. It 
is only by solving the constraint that the Wilson lines 
(and the coset element) appear. This we show below, as 
well as the rewriting in the form of two different metrics 
for (covariant) vector and axial combinations. As an in¬ 
termediate step, we first reinstate the coset elements on 
the IR brane. 

To do this, we go back to the question of symmetry 
breaking by BGs. Rather than directly identifying con¬ 
nections on the IR brane, we use a spurion on the IR 
brane to render the formalism covariant. This amounts 
to performing the identification up to a gauge, as done 
in pI9l| . In the present case, one would introduce a field 
uj (x) living on the IR brane, on which the constraint 
= 0 is applied, lo should in fact be an element of the 
coset of the transformations on the IR brane, i.e. it be¬ 
longs to [SU (tVf) X SU (Nf)] / SU (Nf). In other words, 
it is unitary, and transform as a bi-fundamental under 
the SU (Nf) X SU (Nf) gauge group. We impose the two 


constraints 

0 = (x) 

= (x) - ii?^ (x, li) to (x) + iuj (x) (a;, ^i)(Al) 

0 = L 5 f^{x,li) + uj'^ {x) R 5 fj,{x,li)u; {x), (A2) 

which are both invariant under the 5D SU {Nf) x SU (-/V/) 
gauge group. They provide a covariant replacement for 
the BCs on the IR brane. 

This defines the model as well as BCs: one can 
then choose to gauge away the coset element w, 
the rest of the algebra is then strictly identical to 
that of ^. Alternatively, one may want to pre¬ 
serve covariance all along, in which case the defini¬ 
tions should be slightly modified: for instance, the 
rotations ^l,r(x,z) appearing in the field redefini¬ 
tions used to obtain appropriate vector and axial fields 

Vm,Am = 5 (^M — ITm) ± Cl (^M — ITm) ^l| 

should involve Wilson lines, but also the held lo. The 
constraint to respect is that 

Cfl 2) cl (a:, (x,a;'')^|a;(x) 


idea. We have already mentioned the constraint to be 
imposed 


0 = D^^l{x,z) 

= (x, z) — ii?5 (x, z) fl (x, z) 

-I- in (x,z) L 5 (x,;i), (A4) 

as well as the BC for 


(x, z = h) = 0, (A5) 

when we wrote down the action frm . Let us now de¬ 
scribe how, starting from a lagrangian that respects lo¬ 
cality in the hfth dimension, as well as covariance under 
the 5D SU {Nf) x SU {Nf) gauge group, the constraints 
above induce two different metrics for vector and axial 
held strengths. 

The solution to the constraint (IAtIi can be written as 


n(x, z) = P-|exp^i^ dx®i?5 (x, x^)^ I n (x, U) 


X P 


exp(./Ux%(..x=) 


(A 6 ) 


X 


exp 



dx^L^ (x, X® 


(A3) 


With this dehnition, the matrix U (x) collecting the pions 
is still written as U = (x,U)C^ {x,Iq), and all equa¬ 

tions follow identically as in [J, but now preserving the 
full covariance at each step. 

This shows how a 4D spurion on the IR brane reintro¬ 
duces the coset element necessary for a covariant descrip¬ 
tion of the symmetry-breaking conditions at the bound¬ 
ary. Note that the original lagrangian was local, but our 
held redehnitions are non-local in the hfth dimension, 
since they involve Wilson lines: this will always be the 
case in such a model where the GB mode is the Wilson 
line, and therefore non-local. We now move to the ques¬ 
tion of the bulk auxiliary held Q (x,z): we start from a 
local and covariant lagrangian generalizing the previous 


Using the BC (IA5II then brings us back to the previous 
case, with the replacement lo (x) i—*■ U (x, h), and all the 
equations follow. In particular, we see that 

{wv {z) + WA (2:)) {LmsLnT + RmsRnt) 

-\- 2 {wv — wa) {Lmn^^Rst^) 

= 4wv (VmnVst^ + 4wa (^Amn^st^ , (A7) 

conhrming the interpretion that this term splits the two 
metrics, at the level of the (covariant) vector and ax¬ 
ial combinations of helds. From here on, the extraction 
of the lagrangian in terms of the pions U and the KK 
modes I 4 ,, A„ proceeds as in Q. The pion held is dehned 
again by U (x) = (x, Iq) Cl ^ 0)5 which is also equal 

to Q{x,lo), hence the straightforward identihcation for 
the quark mass terms of Section nTDi 
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APPENDIX B: APPROXIMATE EXPRESSIONS 


We give approximate expressions for some quantities, at linear order in the condensates. Using reduced values of 
the condensates 04 = (O 4 )/(iVc/j"^) and Oq = (Oq) / we have ® 


u = 


VK 




3.83 


y/dirli 

(1 - 3.35o4 + 0.38o6) + O (o^) 


O (o^) 


Ma, 


U.4 — —31^2 — 


- 6 L 1 = - 

Lg = 


(1 + I.OI04 + 0.29o6) + O (o^) 
11 W 


15367r2 

N. 


1457r^ 




Lio — — 


647r2 
Nr 


( 257r^ StT^ \ f 2\ 

(l + ^04 + ^a6j+0(o), 
257r^ 37r^ \ / 2\ 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B 6 ) 


It turns out that the expressions are accurate to a few percent for the range of interest (the one of Figure QJ. 

I- 


Adding fnsii and IIB 6 II shows that, if (Oq) is of the or¬ 
der of magnitude expected from factorization on obtains 
Lg -I- Lio ~ 2 X 10“^, i.e. one order of magnitude be¬ 
low the experimental value. Similarly, the model under¬ 
predicts the decay if factorization is respected. 

This can be understood as follows: a model that repro¬ 
duces the width Tr^aiu^ has Jai ^ 0.13. Assuming that 

® The numerical coefficients stem from integrals of Bessel func- 
tions. 


the sum rule is saturated by the oi one can esti¬ 

mates Va^^TT-y from the knowledge of Lg -|- Lip from fac¬ 
torization. A suppression by an order of magnitude for 
Lg -I- Lio yields a suppression by two orders of magnitude 
for Tai^TT'y- Indeed, we find that ~ fewkeV for 

the favored region of Figure ^ whereas experimentally 
one has raj^ 7 r 7 = 640 ± 246keV 
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